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ABSTRACT 



The magnetic Reynolds number Rm, is defined as the product of a 
characteristic scale and associated flow speed divided by the microphysical 
magnetic diffusivity. For laminar flows, Rm also approximates the ratio of 
advective to dissipative terms in the total magnetic energy equation, but for 
turbulent flows this latter ratio depends on the energy spectra and approaches 
unity in a steady state. To generalize for flows of arbitrary spectra we define 
an effective magnetic dissipation number, RM,e, as the ratio of the advection 
to microphysical dissipation terms in the total magnetic energy equation, 
incorporating the full spectrum of scales, arbitrary magnetic Prandtl numbers, 
and distinct pairs of inner and outer scales for magnetic and kinetic spectra. 
As expected, for a substantial parameter range RM,e ~ 0{1) « Rm- We also 
distinguish i?M,e from i?M,e where the latter is an effective magnetic Reynolds 
number for the mean magnetic field equation when a turbulent diffusivity is 
explicitly imposed as a closure. That RM,e and RM,e approach unity even 
if Rm » 1 highlights that, just as in hydrodynamic turbulence, energy 
dissipation of large scale structures in turbulent flows via a cascade can be 
much faster than the dissipation of large scale structures in laminar flows. This 
illustrates that the rate of energy dissipation by magnetic reconnection is much 
faster in turbulent flows, and much less sensitive to microphysical reconnection 
rates compared to laminar flows. 
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1. Introduction 

Magnetic fields play an observable and dynamical role in a range of astrophysical 
and laboratory plasmas. When the scales of interest are larger than the particle mean 
free path, magnetohydrodynamics provides a suitable approximation (Mofffat 1978, 
Parker 1979, Biskamp 1997). Derived from Maxwell's equations and Ohm's law, the 
induction equation describing the time evolution of the magnetic field B in non-relativistic 
magnetohydrodynamics (MHD) is given by 

a^B = V X (v X B) + AV^B, (1) 

where A is the magnetic diffusivity, and v is the velocity, determined from the fiuid 
momentum equation. For incompressible fiows, the momentum equation incorporating the 
Lorentz force is 

dt\ = -V • Vv - VP/p + (V X B) X B/47rp + i/V^v, (2) 

where P is the pressure and v is the kinematic viscosity, and p is the density. 

Taking the approximate ratio of magnitudes of the second to third terms in (1) 
introduces the magnetic Reynolds number 

Rm,i = Iv/X, (3) 

where I is a macroscopic scale of field variation for the system being studied. When I is the 
largest macroscopic scale of the system, we write Rm, eliminating the subscript I. 

For laminar flows in which Rm » 1, the last term of (1) is often ignored. An 
important consequence of ignoring this term is magnetic flux freezing, / B • dS = 
(Parker 1979) which implies that the integrated product of the magnetic fleld and the 
surface area it penetrates, or the magnetic fiux moving with the plasma, is conserved. 
Dotting (1) with B gives the same Rm for the corresponding terms in the magnetic energy 
equation so hat Rm >> 1 would also imply that the resistive dissipation term is small 
compared to the other dynamical terms. When most of the magnetic energy is on the largest 
scale, the implication is that the dissipation of the total magnetic energy is negligible. 

The above physical implications of Rm >> 1 do not apply if the fiow is turbulent. 
When an initially laminar fiow becomes turbulent, the field and fiow evolve to develop 
non-trivial turbulent spectra (e.g. Biskamp 1997, Biskamp 2003). Eqn. (3), with I taken 
as the largest macroscopic scale of the system, is no longer a good approximation to the 
ratio of the advection to dissipation term in the total magnetic energy equation even if the 
energy is dominated by structures on that scale. The dissipation term is not necessarily 



-3- 



negligible because it depends more strongly on high wavenumber components of the 
turbulent spectrum than do the advection terms. 

The importance of dissipation in MHD turbulence just described, is directly analogous 
to the case of hydrodynamic turbulence. For the latter, the Kolmogorov steady-state 
assumption (Kolmogorov 1941; see also Davidson 2004) that the energy dissipation rate per 
mass is constant on all scales through the cascade implies that the energy dissipation rate 
at the dissipation scales exactly equals that at the input scale. Therefore, the dissipation 
term in the total energy evolution equation exactly balances the nonlinear energy transfer 
terms when all scales are included and a steady-state is assumed. Extracting dimensionless 
measures of hydrodynamic turbulence dissipation rates have confirmed this basic conceptual 
picture (e.g. Pearson et al. 2004). 

In making contact with these concepts from hydrodynamic turbulence, we find it 
particularly important to emphasize that dissipation is non-negligible in MHD turbulence 
even though Rm » 1. This very much contrasts the physical implications commonly 
attributed to laminar Rm » 1 MHD flows. Characterizing the MHD energy depletion 
time scales for large Rm turbulent MHD flows such as those of stellar convection zones, 
accretion disks, and galactic interstellar media, is important for understanding the magnetic 
field evolution, conversion of magnetic energy into particles, and requirements for magnetic 
field sustenance. Since these fiows are characterized by a spectrum of scales (Biskamp 1997, 
Biskamp 2003), the ratio of advection to dissipation terms in the total magnetic energy 
equation is not in general Rm but a generalized magnetic dissipation number RM,e which 
depends on the magnetic and kinetic energy spectra. 

Quantifying RM,e, requires incorporating the spectrum of magnetic and velocity fields 
along with the magnetic Prandtl number Ptm = i^/X- Here we calculate RM,e as a function 
of the magnetic spectral index and Ptm, the latter of which depends on the microphysical 
diffusivities and thus on Rm- Although Rm (as defined with the largest scale) does help 
determine whether the flow would have become turbulent in the flrst place, it is not a 
good indicator of the time scale for overall energy dissipation when the flow is turbulent. 
For regimes in which RM,e is largely independent oi Rm, the total energy depletion rate 
becomes largely independent of the specific small scale mechanism of dissipation; turbulence 
cascades the energy to scales where the dissipation time is short. 

In section 2, we derive the approximations for the corresponding second and third 
terms of (1) in the magnetic energy equation, incorporating the turbulent spectra and 
triple correlations for the case without a mean field. We then plot RM,e as a function of 
the magnetic spectral index and Ptm- In section 3, we distinguish i?M,e from a different 
quantity i?M,e for the mean field equation. In section 4, we conclude with a discussion of 
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the implications of RM,e << 1 and RM,e << Rm- 

2. RM,e in the absence of a mean field 

We write the magnetic field as a sum of mean and fiuctuating contributions such that 
B = B + b where B is the mean field that survives the averaging over scales much larger 
than the scale of any fiuctuations incorporated in b. We use the same approach for the 
velocity, but we will assume here that the mean velocity V = 0. 

We first consider B = for which equation for the total fiuctuating magnetic field 
(written in Alfven speed units), obtained by subtracting the average of (1) from (1), is 

a^b = V X (v X b) - V X (v X b) + AV^b. (4) 

Dotting this equation with b and averaging, we obtain 

-dt{h^) = (b • V X (v X b)) + A(b • V'b). (5) 

The effective magnetic dissipation number RM,e represents an estimate of the ratio of 
magnitudes of the second to third term in (5). 

When an initially laminar system is subject to a source of turbulence in the velocity 
field, both the velocity and magnetic field quickly evolve nonlinearly to acquire turbulent 
spectra. Even if the magnetic initial field were weak, the system typically evolves to a state 
of near equipartition between total kinetic and magnetic energies, but each with distinct 
spectra (Schekochihin et al. 2002, Schekochihin et al. 2004, Haugen et al. 2003, Haugen 
et al. 2004). The triple fluctuation terms in e.g (5) mediate the nonlinear sustenance 
of the turbulent spectra. To obtain quantitative estimates of (10) and (11), we write 
the magnitude of the magnetic field and velocity in terms of angle-integrated kinetic and 
magnetic energy spectra EM{k) and Ev{k) respectively. We have 

l-q 

v(k) = (2kEv(k))'/' = vo (^^) ' (6) 

and 

b{k) ^ {2kEM{k))'/' ^ bJ^] ' , (7) 

where q and m are the kinetic and magnetic energy spectral indices, Vq and bo are the 
magnitude of the velocity and magnetic field on the outer turbulent scale, and /cq,^, /co.b are 
the wave numbers for the outer scales of the assumed kinetic and magnetic spectra. 
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We also allow for different inner scales, or equivalently, allow the magnetic Prandtl 
number Ptm to be a free parameter. To accommodate this, note that at the scale k^^, the 
magnetic spectrum is truncated by dissipation. We assume that the kinetic energy cascade 
is mediated by local interactions so that u ~ Vy/ky, where and k^ are the velocity and 
wave number of the viscous scale. There must also be a characteristic speed v^ff such that 
A = Vf,ff/kx- However, the value v^jf need not be vx (the value of v at the scale k^^)^ and 
cannot be for Ptm >> 1, since then w;^ ~ 0. Larger scale velocities fold the field which can 
show up as magnetic power on small scales (Schekochihin et al. 2002, Schekochihin et al. 
2004, Haugen et al. 2003, Haugen et al. 2004). We then posit that hx ~ v^ff, since hx (the 
field at kx) draws its energy from t'e//- Using the definition of Ptm = X' ^® then have 



X Vy 

A ~ — ~ — — — (8) 
kx Ptm K 



Eq. (8) implies 





RM^^^ = PrM[^\ (P^) =?(f^) . (9) 

Assuming that the magnetic and kinetic spectra for k < Min[kx,ky] are mediated 
primarily by local interactions, we can approximate each nonlinear term on the right of (5) 
by an integral over the magnitude of k. For ko^b > ko^y, the second term of (5) gives 



(10) 



where kd/ko^f, = Min[ky/kQ^b, kx/ko^b] = Min 



k),v ( Rm \ ( RmW) \ '"+1 
h),b KPrMj ' V ''0 



, using (9). The 



inequality in (10) arises because the left side includes the competing effects of turbulent 
amplification and turbulent diffusion of the magnetic field. The norm of the third term of 
(5), similarly averaged, is given by 

ku Jko,b \ko,bJ Rm \ko,vJJi 

where we have again used (9). Because (11) involves only dissipation, whereas the left side 
of (10) could involve amplification or diffusion, RM,e, estimated by the ratio of the right 
side of (10) to the right side of (11), represents an upper limit to the ratio of the left hand 
sides of those equations. 

To use the ratio of (10) to (11) to estimate RM,e: we need values for q, m, bl/v^, Rm 
and PrM- We will assume that q — 5/3 and plot i?M,e as a function of independent variables 
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m and Ptm for various values of Rm- For m > 1, b^/vg can be approximated by the ratio 
of total magnetic to kinetic energy because 

bl _ Em K-'>/^dK ^ Em K-'>/^dK ^ 3Em 1 - (RM/PrM)-'^' Em 
^0 Ev Ev 2Ev -^(i-(R,,hJvn)'^\^ Ev' 

(12) 

where the similarity follows self-consistently for RMbo/vo, RAi/PrM >> 1- For modest Rm 
and Ptm, the same approximation is also satisfactory for m = 1, as the integrals in the 
second term of (12) become ratios of logarithms. In general, 6q is of order Vq for the range 
of m considered. Remember that bo and are defined at A;o,6 and ko^y respectively and the 
latter two wave numbers need not be equal: Although we assume m and q are constant 
over the ranges of < k < k\ and ky^ < k < ki, respectively, the allowance of k^^ ^ kb^ 
and kx ^ k^, provides flexibility in incorporating a range of spectra compatible with MHD 
simulations (Schekochihin et al. 2002, Schekochihin et al. 2004, Haugen et al. 2003, Haugen 
et al. 2004). 

Figs. 1 and 2 show i?M,e calculated from the ratio of (10) to (11) for 10~^ < Ptm < 
and 1 < m < 1.7 using two values of kyfl/kbfl and two values of {bo/voY. For Fig. 1, 
Rm = 10^ and for Fig. 2 Rm = 10^. To ensure k„ » k,,,Q, wc niTist have Rm/Ptm >> 1, 
which explains why we have chosen the upper limits on Ptm to be Rm/^O in each plot. 

The case of {bo/voy — 1/4 used in Fig. Ibd and 2bd is motivated by saturated states of 
non-helical MHD turbulence which show that in cases when the magnetic energy is initially 
weak, it builds up to a fraction of equipartition of the total kinetic energy when the latter 
is steadily driven (Schekochihin ct al. 2004, Haugen et al. 2004, Haugen ct al. 2003, Maron 
et al. 2004, Haugen et al. 2004). Wc have not required a steady state, but in that case the 
two terms on the right of (4) would exactly balance. This could be used to set the value of 
{bo/vo)^ for a given m, ky^o/kb,o, Ptm, and Rm, thereby enforcing RM,e = 1- However, the 
scaling approximation for the first term on the right of (10) is crude, and does not take into 
account that terms within that term can have different signs, reducing its value and that of 
i?M,e- In the steady state, RM,e, could therefore modestly exceed unity using our estimation 
procedure. As seen in Fig. 1 however, the important point is that i?M,e << Rm robustly. 

The dependencies of i?M,e on Rm, Ptm, ky^o/kb^o, and {bo/voY are evident in Figs. 
1 and 2. In each panel, as m increases for a fixed Rm, the magnetic spectrum steepens 
and becomes more heavily weighted by low k for the same value of (bo/vo)^, increasing 
i?M,e- A larger value of {bo/vof' lowers RM,e, as can be seen by comparing the two rows of 
Fig. 1. This is expected: for a given m, increasing {bo/vo)^ increases the upper bound in 
(11) and that integral is more sensitive to the upper bound than the integral of (10). Since 
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RM,e oc {kyfi/kbfi) 2 , increasing /c^,o/^6,o also increases i?M,e as evidenced by the figures. 

Each panel of Fig. 1 shows that for a fixed m, increasing Ptm lowers RM,e- This is 
because more spectral range is available for the magnetic field at large k as Pvm is increased, 
thereby increasing the relative weight of high k contributions of the field in the dissipation 
term. At large Ptm and low m, i?M,e since the shallow spectrum again increases the 
relative importance of high k for the dissipation term; the advection term is ruled by the 
k^ «kx velocity cutoff and is less sensitive to increases in Ptm for Ptm >> 1. 

Comparing Figs. 1 and 2, we see that increasing Rm increases, i?M,e only weakly, 
highhghting the role of the spectra and the universahty of RM,e « Rm- We reiterate that 
i?M,e is calculated from the ratio of the advection to microphysical dissipation terms, NOT 
by approximating any of the advection terms via a turbulent diffusivity as is commonly 
done in mean field formalism (Mofffat 1978, Zeldovich et al. 1983). We discuss the latter in 
the next section. 



3. Eflfective dissipation number for the mean field equation 

When the dynamical equations of a turbulent system such as an accretion disk 
(Shakura and Sunyaev 1973, Frank et al. 1992, Balbus and Hawley 1998) are assumed 
to obey a symmetry such as axisymmetry, these equations describe the evolution of mean 
quantities, not total quantities. Turbulence violates the global symmetry locally and only 
the mean quantities obey the symmetry. Accretion disk theory that invokes a turbulent 
viscosity (Shakura and Sunyaev 1973) in the axisymmetric velocity equation is in fact a 
mean field theory, even when not explicitly presented as such. The same is true for the mean 
magnetic field equation (Mofffat 1978, Zeldovich et al. 1983). The associated turbulent 
diffusivities are formally constructed from closures that replace mean nonlinear correlations 
of macroscopic turbulent fluctuations in the momentum and induction equations by 
diffusion terms, NOT from the replacement of microphysical dissipation terms. This is an 
important distinction from the previous section which involved no closures on the second 
term of (5); there we simply took the ratio of that term to the third term to define RM,e- 

More explicitly, consider the equation for evolution of the mean magnetic field. This 
equation is obtained by averaging (1), which gives (Mofffat 1978, Parker 1979) 

9tB = V X (v X b) + V X (V X B) - V X AVB. (13) 

When the turbulence is assumed to be isotropic, three dimensional, and incompressible, 
and (v X b) is time independent, mean field theory with the minimal r closure (Blackman 
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and Field 2002) gives 

(v X b) ~ roB - r/3V x B + /(V), (14) 

where r ~ is a non-linear turbulent damping time, and in 3-D, t/9 ^ (v^t) acts 

as a diffusion coefficient. The a term, important for the helical dynamo (Mofffat 1978, 
Parker 1979, Brandenburg and Subramanian 2005) survives only when conditions (such as 
rotation + stratification) are able to sustain a mean pseudoscalar. The last term in (14) 
schematically indicates additional terms that are a function of V. 

Here we do not discuss the precise value of a and (5 in (14), rather, we simply want 
to emphasize that these terms come from turbulence and would be absent without it. In 
the simplest case for which the second and fourth terms of (14) vanish but the third term 
remains, we can define the effective magnetic Reynolds number for the mean field equation 
(13) as RM,e = |vx(^^4^)(VxB)| • terms of the velocity spectrum, if A << r/3, we can 
approximate the total diffusion by the latter and bound it by its upper limit such that 

rP + Xr^rP^ (v^t) ~ /'^^'°" K-^'^+'^'dK. (15) 

For q = 5/3, r/5 ~ fo/fco/u- This highlights that in the presence of turbulence, the mean 
field can diffuse at a rate determined by a macroscopic turbulent diffusion coefficient rather 
than the microphysical diffusion coefficient. That is, 

RM,e < ^ = yr^, (16) 

which can be unity for V ^ Vq and Lko^y ~ 1. 

We emphasize again that RM,e of the previous section has a different definition from 
RM,e just defined. The former measures the ratio of advection terms to the microphysical 
dissipation terms in the total energy equation, while the latter applies only to the mean 
field equation and depends on our having introduced a macroscopic turbulent diffusivity 
T/S to replace some triple fiuctuation terms within the advection term by a diffusion 
term. The total mean field diffusion coefficient is sum of turbulent and microphysical 
diffusivities and so the denominator of i?M,e also includes contributions from both. When 
the turbulent diffusivity greatly exceeds the microphysical diffusivity, the the denominator 
of i?M,e is dominated by the former. This contrasts the denominator of i?M,e which only 
involves microphysical dissipation. The reason that i?M,e can also be small, as discussed 
in the previous section, is that small scale structures lead to a substantial microphysical 
dissipation term even if the microphysical diffusion coefficient is small. 

That i?M,e and RM,e can be small and independent of Rm results because turbulence 
cascades large scale structures down to microphysical scales where dissipation is fast. For 
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driven turbulence, the overall rate at which the cascade occurs is fixed by the forcing. As 
long as there is some microphysical mechanism of energy drain at small scales the rate of 
that mechanism does not strongly influence the overall rate of total energy dissipation; the 
system will develop the structures needed to accommodate the cascade rate. 

The extent to which mean held turbulent diffusion in MHD is suppressed below 
the maximum kinematic macroscopic value by the backreaction from fluctuating and 
mean magnetic fields has been a topic of study (Vainshtein and Cattaneo 1992, Cattaneo 
1994, Blackman 1996, , Rogachevskii and Kleeorin 2001, Blackman and Brandenburg 
2002, Brandenburg and Subramanian 2005). Part of the reason for some apparent 
disagreement is that different researchers have made different simplifying assumptions 
(boundary conditions, steady vs. evolving states, 2-D vs. 3-D) and although each may 
present a result consistent with their assumptions the universality of the conclusions still 
needs to be sorted out. The suppression seems to be less dramatic in 3-D. More simulations 
are needed and further detailed discussion is beyond the scope of the present paper. To 
accommodate a range of probabilities we have written an inequality in (16), but the basic 
principles discussed as to why a cascade would dictate i?M,e << 1 and i?M,e < Rm remain. 

As discussed in Sec 1., the principles regarding the importance of the dissipation 
terms derive from, and are analogous to similar circumstances in the purely hydrodynamic 
case, with analogous effective dissipation numbers defined appropriately. In the absence 
of magnetic fields however, there is no backreaction so the hydrodynamic evolution of the 
mean velocity field and its turbulent diffusion are not subject to the issues raised in the 
previous paragraph. 

4. Conclusions 

We have calculated a dimensionless measure of the relative importance of advective vs. 
dissipative terms for the total magnetic energy equation in MHD that applies for arbitrary 
magnetic and kinetic spectra, and magnetic Prandtl number. We have discussed how the 
standard magnetic Reynolds number Rm defined using the largest gradient scale of a system 
does not accurately approximate this ratio for turbulent systems, even if the magnetic 
energy is dominated by large scale structures. Instead, this ratio is best approximated 
by RM,e, a quantity largely independent of Rm for a range of turbulent spectra (though 
reverting to Rm for a laminar flow) and typically satisfles i?M,e ~ 0{1) « Rm- The latter 
results because small scale structures augment the relative importance of the microphysical 
dissipation terms for a turbulent system compared to a laminar system with the same 
microphysical diffusion coefficient. In a steady state, such as those reached when MHD 
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turbulence is steadily forced, RM,e ~ 1- 

We have also discussed how the effective dissipation number RM,e for the mean 
magnetic field equation is defined differently from i?M,e because the mean field equation 
picks out the gradient scale of the mean field only. A low value of RM,e « Rm results not 
from small scale gradients, but because macroscopic turbulent diffusion coefficients exceed 
microscopic diffusivities in closures that approximate nonlinear cascade terms by diffusion. 

While the results that i?M,e << Rm, RM,e « Rm, are not at all surprising when 
interpreting them analogously to steady-state Kolmogorov hydrodynamic turbulence, 
they do highhght the much greater role of dissipation in large Rm turbulent MHD 

flows compared to large Rm laminar flows. In particular, an important implication of 
RM,e « Rm, RM,e « Rm, and the very weak dependence of i?M,e and RM,e on Rm is 
that the rate at which large scale MHD energy is dissipated can be much less sensitive 
to the actual microphysical dissipation mechanism in turbulent flows when compared to 
laminar flows. On a macroscopic dynamical time, the MHD energy cascades to scales small 
enough such that it can be drained at multiple sites. If the turbulence is steadily forced, 
then the system adjusts to produce enough small scale structures to accommodate the 
driving energy input rate. The overall energy dissipation rate can be insensitive to the 
microphysical reconnection (or other microphysical energy conversion rate) at any one site. 
Instead, the aggregate of dissipation sites emerges to drain the total magnetic energy on 
macroscopic dynamical time scales. Since many astrophysical settings are MHD turbulent, 
this highlights the importance of being precise when asking and studying the question "Is 
magnetic reconnection fast in astrophysics?" 

The principle that small scale structures increase the effective reconnection rate is 
consistent with rapid reconnection mechanisms in the MHD models of Hendrix et al. (1996), 
Galasgaard & Nordlund (1996), Lazarian & Vishniac (1999), and in the rapid reconnection 
phase of Fan et al. (2004). Reconnection studies that start with a large scale laminar field 
and reconnect on coUisionless scales (e.g. Birn et al. 2001, Shay et al. 2007) address a 
different regime. The latter can however apply to each of the small scale multiple dissipation 
sites within a globally turbulent MHD system. In some systems, coUisionless individual 
structures are observationally resolved (e.g. the Earth's magnetosphere) and studying the 
plasma instabilities and plasma turbulence for different conditions on such scales is needed 
to understand when reconnecting structures incur a reconnection rate that depends only on 
its macroscopic properties (Kulsrud 2001). That being said, if a stellar or accretion disk 
corona is considered as a global entity, its overall energy dissipation rate in a steady state 
is determined by the dynamical rate of magnetohydro dynamic energy injection from below; 
the small scale structures develop to accommodate this rate. The overall dissipation rate 
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ls best characterized by Rm « Rm and RM,e << Rm even if the rate of each individual 
dissipation site were to depend on Rm- 
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(bo/Vo)'=l/4,kv,o/kb,o = l/5 
(bo/Vo)'=l/100,kv,o/kb,o = l/5 /P^-^. 




Figure 1: RM,e vs. magnetic spectral index m and magnetic Prandtl number Ptm for 
Rm = 10^. a) {bo/voy = 1/100, K,o/h,o = 1/5 b) {bo/vof = 1/4, K,o/h,o = 1/5 c) 
(^o/'^o)^ = 1/100, K^o/khfi = 1 d) (&o/'^o)^ = 1/4, K,o/kb,o = 1 Note the weak dependence of 
RM,e on 
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(bo/Vo)^=l/100,kv,o/kb,o-l/5 
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Figure 2: Same as Fig. 1 but with Rm = 10^ 
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